AD  742275 


280 


SOME  APPLICATIONS 
OF  A  MIXED  SIGNAL  PROCESSOR 


R.  H.  S  HUM  WAY 

CONSULTANT  TO  THE  SEISMIC  OATA  LABORATORY 


B  JANUARY  1872 


'  hftni  hr 

AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 

Wnkinftm,  D.C. 

f 

Pro/nt  VELA  UNI ,  JRM 


AOVANCEO  RESEARCH  PROJECTS  ABENCY 
Nvclttr  Mmitormf  Rttitrch  Offiet 
ARPA  Oritr  Nt.  1714 


^TTELEDYNE  GEOTECH 

Reproduced  by  ALEXANDRIA  LABORATORIES 

NATIONAL  TECHNICAL 
fNFORMATION  SERVICE 

Springfield,  Va.  22151 

trnmi  f»» fihic  ihmii;  iistiiietim  uu«ite»%  . ) 


*****  hwt  A«to(a  4|M*r  H'  (to  *"  *«•  Tmlmiml  Cmt*  will  to  raa*aaatoto  tor  totomatoa 

mfr'rr1  tonto  wkkk  to*  toM  aawtoto  to  *tor  mitotoi  ar  aaarraarara.  to  (Mi  toMM*H<to>  *•  tow  rwtotoa  aa  Mr  to 
itrtrrrrrr  Jto  rtoi  to  wtotowtom  ftoaa  to  ito  aartora  to  atoato  to  to  toWtoMto  *•  mmmrAr  totoMMrtof  »to 
•tototo/  a*toto».  f/ftor  in/mrr"  ar  topMto,  to  (to  «*aaato  toaaarti  AnjtUt  toto  (to  A/r  faraa  farto/aa/  Jaatoa'toaa  Caaiar. 
ar  (to  H  S  liwwawl. 


t 


Security  Classification 


■“  DOCUMENT  CONTROL  DATA  -  RAD 

(Saeurity  datamation  ol  titl a.  body  ol  abmtract  and  inditing  annotation  mull  bo  antarad  whin  thi  ovirill  riport  is  elaiilliad) 

1  ORIOINATIN  <5  ACTIVITY  (Corponf  Author.) 

TELEDYNE  GEOTECH 

2  •  nCPORT  SECURITY  CLASSIFICATION' 

Unclassified 

ALEXANDRIA,  VIRGINIA 

2  t>  GROUP 

1  REPORT  TITLE 

SOME  APPLICATIONS  OF  A 

.  MIXED  SIGNAL  PROCESSOR 

4  DESCRIPTIVE  NOTES  (Typo  ol  riport  and  Inelualva  data) 

Scientific 

S  AUTHOR)  (Lad  nama,  Itnl  nima.  Initial) 

Shumway,  R.H. 

Consultant  to  the  Seismic  Data  Laboratory 

•  REPORT  OATE 

5  January  1972 

7  •  TOTAL  NO  OF  PAGES 

48 

7b.  NO  OF  REFS 

16 

|  a  CONTRACT  OR  GRANT  NO. 

F33657-72-C-0009 

6  PROJECT  NO 

VELA  T/ 2706 

ta  ORIGINATOR'S  REPORT  NUMBCRfSJ 

280 

ARPA  Order  No.  1714 
<t  ARPA  Program  Code  No.  2F-10 

tb  othic.  report  MOfl)  (Any  othir  numbara  that  may  bi  aaalfnad 
thla  riport) 

10  AVAIL  ASIL  IT  Y /LIMITATION  NOTICES 

APPROVED  FOR  PUBLIC 

RELEASE;  DISTRIBUTION  mWITEI, 

II  SUPPLEMENTARY  NOTES 

12  SPONSORING  MILITARY  ACTIVITY 

Advanced  Research  Projects  Agency 
Nuclear  Monitoring  Research  Office 
Washington.  D.  C. _ 

II  ABSTRACT 


A  technique  to  process  array  data  has  been  developed  to 
determine  when  a  second  signal  is  hidden  in  the  coda  of  a 
primary  signal;  and  to  estimate  the  waveform  of  the  two  signals. 
The  effectiveness  of  the  mixed-signal  processor  has  been  demon¬ 
strated  by  operating  on  various  possible  mixed  signals  formed 
from  recordings  of  earthquakes  at  Tonga  in  the  South  Pacific 
and  at  the  Fox  Islands  in  the  Aleutians.  For  small  arrays  the 
processor  is  found  to  be  substantially  superior  to  simple 
beamforming.  The  distortion  introduced  by  the  use  of  finite 
time  filters  is  found  to  be  negligible. 


14  Kt  V  WORDS 

Signal  Estimation 
Maximum  Likelihood 


Unclassified _ 

Security  Classification 


SOME  APPLICATIONS  OF  A  MIXED  SIGNAL  PROCESSOR 
SEISMIC  DATA  LABORATORY  REPORT  NO.  280 


AFTAC  Project  No. : 

Project  Title: 

ARPA  Order  No.: 

ARPA  Program  Code  No.: 

Name  of  Contractor: 

Contract  No.: 

Date  of  Contract: 

Amount  of  Contract: 
Contract  Expiration  Date: 
Project  Manager: 


VELA  T/2706 

Seismic  Data  Laboratory 
1714 
2^-10 

TELEDYNE  GEOTECH 

F33657-72-C-0009 

01  July  1971 

$  1,314,000 

30  June  1972 

Royal  A.  Hartenberger 
(703)  836-7647 


p.  0.  Box  334,  Alexandria ,  Virginia 


APPROVED 


for  public  RElt*»: 


DISTRIBUTOR  URLIDtTEB. 


ABSTRACT 


A  technique  to  process  array  data  has  been  developed 
to  determine  when  a  second  signal  is  hidden  in  the  coda 
of  a  primary  signal  and  to  estimate  the  waveforms  of  the 
two  signals.  The  effectiveness  of  the  mixed  signal  pro¬ 
cessor  has  been  demonstrated  by  operating  on  various 
possible  mixed  signals  formed  from  recordings  of  earth¬ 
quakes  at  Tonga  in  the  South  Pacific  and  at  the  Fox 
Islands  in  the  Aleutians.  For  small  arrays  the  processor 
is  found  to  be  substantially  superior  to  simple  beam¬ 
forming.  The  distortion  introduced  by  the  use  of  finite 
time  filters  is  found  to  be  negligible. 


TABLE  OF  CONTENTS 


Page  No. 

ABSTRACT 

INTRODUCTION  x 

GENERAL  THEORY  6 

WINDOW  THEORY  FOR  TWO-SIGNAL  FILTERS  9 

ESTIMATION  AND  DETECTION  CAPABILITIES  AT  TFO  13 

SUMMARY  18 

REFERENCES  1Q 


LIST  OF  FIGURES 


Figure  Title 

Array  configuration  for  TFO 


Figure  No. 
1 


Comparison  of  distortion  introduced  by  trun¬ 
cated  optimum  filters  and  ordinary  beamforming. 

19  Channels  at  TFO  with  Tonga  as  S^C.)  and  Fox 
as  S 2  ( . )  . 

Comparison  of  distortion  introduced  by  trun¬ 
cated  optimum  filters  and  ordinary  beamforming. 
Seven  channels  (Z1-Z7)  with  S,(.)  at  17  km/sec, 

308  degrees  and  S 2 ( • )  at  23  km/sec,  300  degrees. 

Comparison  of  distortion  introduced  by  trun¬ 
cated  optimum  filters  and  ordinary  beamforming. 
Seven  channels  (Z1-Z7)  with  S..(.)  at  17  km/sec, 

308  degrees  and  S 2 C • )  23  kin/ sec,  240  degrees. 

Comparison  of  distortion  introduced  by  trun¬ 
cated  optimum  filters  and  ordinary  beamforming  . 
Three  channels  at  TFO  with  Si(.)  20  km/sec, 

0  degrees  and  S2(0  at  10  km/sec,  20  degrees. 

Comparison  of  distortion  introduced  by  trun¬ 
cated  optimum  filters  and  ordinary  beamforming. 
Four  channels  at  TFO  with  Si(.)  at  20  km/sec, 

0  degrees  and  S 2 C • )  at  10  km/sec,  30  degrees. 

Analysis  of  Tonga  Si(t)  event  using  N  =  19 
channels  recorded  at  TFO  at  high  signal  to 
noise  ratio. 

Analysis  of  Tonga  Si(t)  -  Fox  S2(t)  mixture 
using  N  =  19  channels  recorded  at  TFO  at  a 
high  signal  to  noise  ratio. 

Analysis  of  Tonga  Sii(t)  and  Fox  S2(t)  noisy 
mixture  using  N  =  19  channels  recorded  at  TFO. 

Analysis  of  Tonga  S^t)  and  Fox  S2(t)  mixture  for 
N  =  19  channels  showing  bias  in  beam  formed 
estimate  for  Fox. 

Results  of  a  least-squares  technique 
for  estimation  of  one  signal  in  the  pre¬ 
sence  of  another. 


2 


3 


4 


5 


6 


7 


8 

9 


10a 


10b 


LIST  OF  FIGURES  (Cont'd.) 


Figure  Title  Figure  No. 


11 


12 

Maximum  likelihood  estimate  for  Tonga  using 

Tonga-Fox  mixture  in  N  =  7  channels  at  TFO 

(Figure  8)  with  various  possible  azimuths 

for  Tonga.  13 

F  statistics  for  detecting  Tonga  in  Tonga-Fox 

mixture  (N  =  7  center  elements  at  TFO) .  14 

Maximum  likelihood  estimate  for  Tonga  using 

Tonga-Fox  mixture  in  N  *  7  noisy  channels  at 

TFO  with  various  possible  azimuths  for  Tonga.  15 

F  statistic  for  detecting  Tonga  in  Tonga-Fox 

noisy  mixture  (N  =  7  center  elements  at  TFO).  16 


Frequency  wave  number  analysis.  N  =  7  channels 
measuring  Tonga  and  Fox  at  TFO  (1.25  Hz). 

Frequency  wave  number  analysis.  N  =  7  noisy 
channels  measuring  Tonga  and  Fox  at  TFO 
(1.25  Hz). 


LIST  OF  TABLES 


Table  Title  Table  No. 

Values  of  F  Testing  for  the  Presence  of 

Fox  in  Data  Containing  Only  Tonga  1 

Values  of  F  Testing  for  the  Presence  of 

the  Fox  Island  Signal  in  Data  Containing 

Both  Fox  Island  and  Tonga  Signals.  II 

Values  of  F  Testing  for  the  Presence  of 

the  Fox  Island  Signal  in  Noisy  Data  Containing 

Both  Fox  Island  and  Tonga  Signals.  HI 

Values  of  F  Testing  for  the  Presence  of 

the  Tonga  Signal  in  Data  Containing  Both 

Tonga  and  Fox  Island  Signals.  IV 


INTRODUCTION 


A  number  of  complications  are  introduced  into 
signal  estimation  and  detection  procedures  when  one  is 
willing  to  admit  the  possibility  of  interfering  wave¬ 
forms.  This  may  occur  when  a  propagating  noise,  genera¬ 
ted  perhaps  by  a  storm,  appears  on  a  seismic  record 
simultaneously  with  a  propagating  signal  from  an  earth¬ 
quake  or  explosion.  It  may  occur  when  signals  from  an 
explosion  (accidently  or  purposefully)  coincide  with 
those  from  an  earthquake  of  an  appropriate  magnitude, 
depth  and  location.  Either  of  the  above  phenomenon  causes 
one  to  re-evaluate  the  applicability  of  the  single 
signal  plus  noise  model 

Hi:  Y.(t)  =  Sx(t  -  Tj  x)  +  njCt) 

for  j  =  1,  ...,N  time  series  sampled  at  T  points 

t  =  0,  1 . T-l,  where  the  plane  wave  signal  suffers 

t  h 

a  time  delay  T^  at  the  j  sensor.  In  (1)  when  the 
noise  processes  are  independent  (i.e.  n^  (t)  and  n^(t) 
are  independent  for  j  f  k)  with  equal  autocorrelation 
functions,  it  is  well-known  (Kelly,  1965),  (Shumway  and 
Dean,  1968)  that  a  simple  beam  produces  the  minimum 
variance  linear  undistorted  estimate  of  the  signal. 
Futhermore,  an  approximation  to  the  likelihood  detec¬ 
tor  can  be  derived  under  either  of  two  assumptions, 
namely 

(a)  The  noise  series  are  stationary  Gaussian  and 
their  noise  autocorrelation  (spectrum)  is  known  exactly. 
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(b)  The  noise  series  are  stationary  Gaussian  and 
their  noise  autocorrelation  (spectrum)  is  unknown. 

In  case  (a)  the  detection  statistic  is  proportional 
to  the  ratio  of  the  beam  power  to  the  known  noise  power 
and  has  asymptotically  a  chi-square  distribution  with 
2BT  degrees  of  freedom  where  B  is  the  bandwidth  and  T 
is  the  sampling  interval.  If  the  data  are  assumed  to  be 
stationary  over  a  long  period  before  the  signal  arrives, 
the  known  noise  spectrum  may  be  replaced  with  an  esti¬ 
mated  one  and  the  resulting  test  statistic  will  con¬ 
verge  in  distribution  to  a  chi-square  variable  with  2BT 
degrees  of  freedom  as  before.  This  procedure  is  basic 
to  the  on-line  detector  in  operation  at  LASA  (Kobayashi  and 
Welch,  1970).  It  should  be  noted  that  the  consequences 
of  a  change  in  the  noise  spectrum  within  the  signal 
window,  e.g.  a  seismometer  noise  burst,  may  be  serious 
if  this  detector  is  employed. 

If  assumption  (b)  is  made,  one  treats  the  spectrum 
of  the  noise  series  as  an  unknown  parameter  in  the 
likelihood  equations  so  that  an  estimate  of  the  noise 
spectrum  is  made  within  the  signal  window.  One  finds 
this  procedure,  in  the  time  domain,  described  in  early 
works  of  Melton,  et  al.  (1957)  and  Booker  (1965). The 
time  domain  representation  necessitated  the  assumption 
that  n.(t)  was  a  white  noise  process.  Failure  of  this 
method3 initially  can  now  be  ascribed  to  a 
combination  of  instrumentation,  small  arrays  resulting 
in  correlated  noise,  and  non-white  noise  over  the  band 
of  interest.  The  method  is  re- formulated  in  the  fre¬ 
quency  domain  by  Shumway  (1970),  (1971)  and  Shumway 
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and  Husted  (1970)  where  it  is  shown  that  the  ratio  of 
the  beam  power  to  the  estimated  noise  spectrum  converges 
in  distribution  to  a  non-central  F  distribution  with  2B1 
and  2BT  (N-l)  degrees  of  freedom  and  a  non-centrality 
parameter  proportional  to  the  signal  to  noise  ratio  and 
number  of  time  series,  N.  This  enables  one  to  approxi¬ 
mate  the  signal  detection  probabilities  for  a  fixed  false 
alarm  rate  as  a  function  of  signal  to  noise  ratio.  The 
F  detector,  coined  by  Melton,  the  Fisher  detector,  has 
been  applied  by  Blandford  (1970) at  TFO;  he  provides 
computational  procedures  for  setting  the  theoretical 
detection  performance  characteristics  and  shows  that 
results  can  be  achieved  which  compare  favorably  with 
those  of  skilled  analysts  reading  film  data.  Cases 
where  off  beam  events  fail  to  give  false  alarms  with 
the  Fisher  detector  are  shown. 


An  extension  of  the  model  implied  by  assumption 
(a)  was  given  by  Capon  et  al.  (1967)  who  assumed  that 
the  noise  was  correlated  between  sensors  implying  a 
known  spectral  matrix.  The  approximate  likelihood  esti¬ 
mate  for  the  signal  in  this  case  is  weighted  by  con¬ 
volution  with  a  vector  impulse  response  function  pro¬ 
portional  to  the  inverse  of  the  known  spectral  matrix. 
The  potential  effectiveness  of  this  method  is  blunted 
somewhat  by  the  lack  of  a  correlation  in  a  well  spaced 
array  and  the  difficulty  of  forming  an  estimate  of  the 
spectral  matrix  which  does  not  change  over  time  for  a 
closely  spaced  array. 

An  example  of  a  multiple  signal  model  is  the  two 
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(2) 


signal  model 

H2:  Y.(t)  =  SjU-T^)  +  S2(t-T.2)  +  nj(t) 

with  S„(t)  a  second  propagating  plane  wave.  Develop¬ 
ment  of  multiple  signal  models  proceeded  somewhat 
more  slowly  due  in  part  to  the  lack  of  a  fast  Fourier 
transform  algorithm  and  in  part  to  a  failure  to  realize 
that  most  multiple  signal  models  can  be  treated  as 
variants  of  the  time  series  regression  model  of 
Bendat  and  Piersol  ,(1966) .  In  the  multivariate  generali¬ 
zation  (Shumway  and  Dean,  1968)  of  this  model,  we 
regard  (Y .  (t)  ,  j  =  1,...N)  as  a  collection  of  output 
series  related  to  a  collection  of  input  series  (deter¬ 
ministic  functions)  through  a  collection  of  impulse 
response  functions  (signals)  (S^  (t)  ,  j  =  1,...,P)  which 
are  to  be  estimated.  The  appropriate  model  is  for 
(t  =  0,  +_  1,  2...) 

03  P 

Y.  (t)  =  l  I  x.k(t-u)Sk(u)  +  n .  (t) 

J  u=-°°  k=  1  J  J 

where  the  choice  P  =  2,  Xjk(t)  =  6(t-T^k) 

j  =  1 . N,  K  =  1,  2  (6(t)  =  1 

for  t  =  0  and  zero  otherwise)  reduces  (3)  to  equation  (2). 
An  early  application  of  this  technique  by  Dean  (1966) 
was  later  described  in  Shumway  and  Dean  (1968) .  Simi lar 
solutions  to  the  estimation  problem  posed  by  the  special 
case  of  (3),  which  yields  the  two-signal  model  (2), 
can  be  found  in  Schweppe  (1968)  and  Kobayashi 
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and  Welch  (1970).  The  asymptotic  detection  theory  for  the 
case  of  P  signals  present,  according  to  the  signal  and 
noise  model  (3),  is  given  in  Shumway  (1970)  where  it  is 
shown  that  the  ratio  of  signal  power  to  noise  power 
appropriate  for  testing  hypothesis  llj  given  by  (1) 
against  H,  given  by  (2),  converges  in  distribution  to 
a  non-central  F  distribution  with  2BT  P  and  2BT(N- 
degrees  of  freedom. 

In  this  report,  we  will  analyze  in  more  detail  the 
detection  and  estimation  capabilities  of  the  two-signu 
model.  This  will  include  a  case  where  two  signals 
arriving  at  TFO  are  mixed  at  high  and  low  signal  to 
noise  ratios.  A  large  array  consisting  of  19  elements 
and  two  small  arrays  containing  seven  and  three  elements 
respectively  are  considered  for  signals  from  Fox  Island 
in  the  Aleutians  and  from  Tonga  Island  in  the  South 
Pacific.  We  investigate  the  distortion  introduced  using 
the  finite  time  truncated  filters  by  calculating  and 
displaying  the  "window"  through  which  the  true  signals 
are  viewed.  Various  velocities,  azimuths  and  subarray 
choices  are  evaluated  both  with  respect  to  the  distor¬ 
tion  introduced  and  with  respect  to  the  detection  per¬ 
formance.  An  analysis  of  the  two  events  is  made  with 
the  second  signal  assumed  to  be  present  at  various 
fixed  azimuths  and  velocities.  This  indicates  that  the 
two-signal  model  could  be  used  as  a  possible  substitute 
for  a  frequency  wave  number  plot  if  a  primary  signal 
noise  source  of  fixed  azimuth  and  velocity  could  e 
identified  as  the  first  signal. 
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GENERAL  THEORY 


To  test,  the  hypothesis  that  the  second  signal  is 
absent  in  the  model  given  by  equation  (2)  i.e.  (2) 
against  (1),  we  need  the  likelihood  estimates  for  S^t) 
and  S2(t)  under  and  H^.  Consider  first  the  solution 
under  H2,  say 


++  N  T-l 

si  -  ,1  l  hik(t-u)Y  (u)  (4) 

J  k=l  a=0  JK  K  J 

where  j  =  1,2.  Suppose  that  the  zero  mean  noise  process 
nj (t)  is  weakly  stationary  and  Gaussian  with  a 
continuous  bounded  spectrum  determined  by 

TT 

pnnM'i“(t't')  K  (5) 


where  l’nn(.)  denotes  the  power  spectrum  of  the  noise 

assumed  to  be  the  same  for  each  j  =  1 . N.  Furthermore, 

nk(t)  is  assumed  to  be  uncorrelated  with  n.(t).  Now, 
using  the  fact  that  the  two  signal  model  (5)  is  a  special 
case  of  (3),  the  transformed  version  of  the  filter  in  (4) 
is  given  (Shumway  and  Dean,  1968)  by 

Hlk(«)  =  A-1(u))  (NelwTkl  -  A(u) e lwTk2) 

H2k(w)  =  A  1(w)  (NeiwTk2  -  A*  (w)ciwTkl) 

(6) 
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with 


A(w)  *  -  |  A  (to)  | 


(7) 


and 


A(w) 


N 

l 

j-1 


l  eiu<Tjl  •  V 


Kc  take  U..(0)  -  ll2k(0)  -  (2N)"1  to  eliminate  the  singu- 
larity  introduced  by  A(0>  -  0.  Under  hypothesis  II,  the 
estimates  arc  S 2 ( ^ )  *  0  ani* 


(8) 


S^ct)  ■  N*1  j  Y,(t  .  T.,) 
1  j  "  1 


(9) 


The  form  of  the  likelihood  dctccvor  for  testing  II, 
against  II,  depends  upon  (Shumway,  1970)  reconstructing 
the  noisc“traccs  at  each  level  using  the  estimated 
signals  under  II,  and  II,.  This  implies  that 


n*(t)  ■  Yj(t)  -  S^(t-Tj,) 


(10) 


and 


"j 


(t) 


Y.(t) 


-  sr^-v 


(t-Tj2) 


If  we  then  define  the  estimated  spectra  of  the  two 

predicted  noise  processes  as  Pnjnj and  i  njnj(“,) 
respectively,  the  average  noise  spectra  under  II,  an 


(ID 
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II ->  would  he 


N 


'nn(“>  *  S‘'  jl,  l'nf  n}  <“> 


and 


'■.»>>  •  *■'  j,  •■;*.(-) 

j*1  j  j 


(12) 


(13) 


Suppose,  then,  that  a  subset  of  I.  frequencies  (a  band 
of  width  B)  about  the  point  w  can  be  found  such  that 

^'nn1  *'  nni*  *'nn^#^  a*’c  approx imately  constant  over  that 
subset  (band).  Then,  an  l:  statistic  with  2BT  and  2BT(N-2) 
degrees  of  freedom  is  given  approximately  by 


!’„(»)  -  V  (w) 

l:(2»T,2BT(N-2))  5  — - Oil _  (K-2) 

P«W 

1 1  1 1 1  is  true  where  l’nn(*>)  and  l’nn(f*»)  now  represent 
spectral  estimates  smoothed  over  the  bandwidth  of 
interest,  the  numerator  of  (14)  measures  the  improve¬ 
ment  in  going  from  the  one  signal  to  the  two  signal 
model.  If  the  alternate  hypothesis  II,  is  true,  (14) 
becomes  a  non-central  I-  with  non-centrality  parameter 
proportional  to  A(w)  and  the  average  noise  power  spec- 
trim  of  the  second  signal  over  the  band  of  interest. 

Before  proceeding  to  detailed  examples,  involving 
seismic  data,  we  examine  the  resolution  and  bias  of  the 
truncated  time  maximum  likelihood  filters. 


(14) 


-8- 


WINDOW  THEORY  FOR  TWO- SIGNAL  FILTERS 

Since  the  Maximum  likelihood  filters  (5)  are  un¬ 
biased  only  for  the  case  where  infinite  two-sided 
operators  can  be  constructed,  it  is  reasonable  to  ask 
how  well  the  truncated  time  versions  reproduce  the 
signal  of  interest.  This  was  accomplished  by  noting 
that  the  expected  value  of  any  filter  output  which 
estimates  S^(t)  (see  (4)  for  example)  may  be  expressed 
as 

T- 1  T- 1 

I  (§ .  (t ) )  ■  l  !:n(t-u)S.(u)  ♦  l  F,.(  t-u)S2(u)  (IS) 

J  u-o  J  u»0 

for  j  •  1,  2.  For  example,  with  j  ■  l.Fjj(t)  is  the 
transfer  function  which  relates  the  true  signal  1  to 
the  estimated  signal  1  and  l^jft)  is  the  transfer 
function  showing  the  amount  of  signal  2  which  leaks 
into  the  estimate  for  signal  1.  In  the  ideal  case 
Fjj(t)  would  be  a  unit  spike  at  zero  and  l:2j(t)  would 
be  zero  for  all  t.  The  deviation  from  these  idealized 
values  for  the  "window"  functions  Fjj(.)  and  l:2j(.) 
measures  the  distortion  or  bias  of  the  filters. 

Similarly,  F12(t)  is  the  amount  of  signal  1  leaking 
into  the  signal  2  estimate  and  F-^CO  the  wwownt  ol 
signal  2  appearing  in  the  estimate  for  signal  2.  The 
2x2  matrix  of  time  functions  characterizes  the  bias 
of  the  filters  in  the  same  way  that  the  spectral 
window  function  characterizes  the  bias  of  the  various 
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methods  for  estimating  the  power  spectrum. 

in  order  to  compare  the  bins  of  the  maximum  likeli¬ 
hood  filters  with  the  bins  of  ordinary  bens,  forming  over 
n  reasonable  range  of  initial  conditions,  an  array  ol 
19  elements  at  TFO  (Figure  1),  was  chosen  and  it  was 
assumed  that  signals  were  arriving  with  vc'ocities 
near  20  km  per  second  on  atim.ths  distr.butcd  about 

the  circle. 

As  a  first  test  ease,  it  was  assumed,  as  in  sub- 
sequent  examples  to  be  given  later,  that  signals  from 
Fox  Island  in  the  Aleutians  and  Tonga  in  the  South 
Pacific  arrive  at  all  10  elements  and  arc  f.ltered 
under  !!,(•)  (ordinary  beam  forming)  and  llj(**)  '.maximum 
likelihood).  Figure  2  shows  the  disto  • ion  introduced 
and  we  note  that  the  maximum  ikelihood  procedure 
produces  virtually  no  distortion.  In  the  beam  forming, 
the  two  signals  are  not  distorted  by  their  own  wave¬ 
forms  but  signal  1  receives  a  small  component  from 
signal  2  and  vice  versa.  One  could  envision  severe 
distortion  only  in  the  ease  where  the  amplitude 
of  signal  2  is  high  relative  to  signal  1  or  vice  versa. 
The  signals  a. c  well  separated  in  velocity  and  acimuth 
so  that  one  would  not  expect  problems  using  either 

procedure . 

in  order  to  examine  distortion  on  a  somwhat  smaller 

arrav,  the  seven  center  elements  “  =7  in  1  igurc  l’ 

were  chosen  and  it  was  assumed  that  signal  1  arrive  at 
,  veiocity  of  17  km/scc  and  azimuth  of  308  degrees.  The 
second  signal  was  assumed  to  have  a  velocity  of  23  km/aec 
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figure  1.  Array  configuration  for  T10. 


and  a  collection  of  19  different  azimuths  running  from 
0  to  360  degrees  were  investigated.  In  general,  the 
results  were  the  same  as  for  Figure  2  except  that  the 
side  lobes  of  the  beam  distortion  are  1/8  the  peak.  Two 
interesting  worst  eases  appear.  The  first  ease  in 
Figure  3  shows  the  distortion  when  the  second  signal  is 
at  300  degree*  azimuth,  i.c.  separated  8  degrees  in 
azi.iuth  and  6  km/scc  in  velocity  from  signal  1.  The 
maximum  likelihood  filters  arc  relatively  undistorted 
in  this  ease  but  the  bcamforming  filters  show  a  definite 
bias.  A  ease  where  the  likelihood  filters  give  more  bias 
is  shown  in  Figure  4  although  one  can  see  that  the  bias 
introduced  is  still  negligible  when  compared  with  that 
introduced  by  the  beam.  Figure  4  corresponds  to  the 
approximate  velocities  and  azimuths  of  the  Tonga  and 
Fox  events. 

In  order  to  see  how  far  it  is  possible  to  go  before 
the  distortion  gets  large,  the  three  center  elements  zl, 
z2  and  z3  were  investigated  for  a  number  of  different 
velocities  and  azimuths.  Figure  5  shows  a  ease  where  the 
two  signals  arc  separated  by  10  km  per  second  in  velocity 
and  30  degrees  in  azimuth.  We  see  that  the  likelihood 
filters  still  produce  a  relatively  undistorted  signal. 
Figure  6  shows  the  worst  distortion  ever  produced  by  the 
likelihood  filters.  In  this  ease,  the  elements  z8,  zlO, 
zl4  and  zl6  were  investigated  at  the  sane  velocity  azimuth 
separation  as  in  Figure  4,  It  is  evident  that  some  leak¬ 
age  can  occur  in  the  worst  eases. 

To  summarize  the  results  of  the  preceding  discussion, 
we  remark  that  for  most  combinations  of  velocities  and 
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figure  4.  Comparison  of  distortion  introduced  by  truncated 
optimum  filters  and  ordinary  heamforming.  Seven  channels 
(Z1-Z7)  with  Sj(.)  at  17  km/sec,  308  degrees  and  at 

23  km/sec,  240  degrees. 


II  d 


sf  ft)  IN  s/+ft) 


sf  ( t )  IN  St++W 

OPTIMUM 

S2  It)  IN  s**ft)  S2  ft)  IN  S2++(t) 


Sf  ft)  IN  S*  ft)  Sj  ft)  IN  S+(t) 

BEAM  FORMED  i 


S2  ft)  IN  S*ft)  S2  ft)  IN  S2+(t) 


Figure  5.  Comparison  of  distortion  introduced  by  truncated 
optimum  filters  and  ordinary  beamforming.  Three  channels  at 
TFO  with  S 1 C . 3  tO  km/scc,  0  degrees  and  Si(.)  at  10  km/sec, 
20  degrees. 
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azimuths  the  distortion  introduced  by  the  maximum 
likelihood  filters  is  remarkably  small.  The  limit  on 
the  techniques  performance  in  practice  will  probably 
arrive  from  lack  of  signal  correlation  from  sensor  to 
sensor.  This  might  cause  the  contamination  between 
signals  to  be  greater  than  that  shown  in  Figures  2-6. 

The  window  computation  of  this  section  provides  a  method 
for  evaluating  the  signal  separating  capability  of  a 
subarray  relative  to  two  plane  waves  propagating  with 
a  given  velocity  and  azimuth.  In  fact,  if  certain  sub¬ 
arrays  are  chosen  for  monitoring  areas  which  might  be 
conducive  to  hiding  explosions  in  earthquakes,  the 
window  computation  can  determine  the  proper  location, 
configuration  and  spacing  of  the  subarray  needed  to 
resolve  the  two  signals.  An  extension  of  the  technique 
to  take  account  of  the  loss  of  signal  correlation  with 
distance  would  be  needed  to  design  the  optimum  array  in 

practice. 
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ESTIMATION  AND  DETECTION  CAPABILITIES  AT  TFO 


In  this  section  we  investigate  in  more  detail  a 
test  set  of  data  consisting  of  a  mixture  of  an  event 
I rom  Fox  Island  and  an  event  from  Tonga,  The  real 
signals  are  mixed  with  real  noise  from  TFO  at  high  and 
low  signal  tc  noise  ratios.  The  velocities  and  azimuths 
arc  approximately  17  km/sec  and  308  degrees  for  Fox 
Island  and  23  km/scc  and  241  degrees  for  Tonga,  It  is 
frequently  convenient  to  interchange  the  numbering 
of  the  Tonga  and  Fox  signals,  so  a  certain  amount  of 
caution  should  be  exercised  in  reading  the  following 
discussion. 

Consider  first,  a  case  where  19  channels  of  data 
are  available  at  a  very  high  signal  to  noise  ratio. 
Figure  7  shows  a  case  where  only  the  Tonga  Island 
(in  this  case  ( . ) )  is  present  at  a  high  signal  to 
noise  ratio.  The  19th  channel  is  shown  along  with  the 
maximum  likelihood  estimates  Sj4(.)  and  S^C.)  under 
**2*  Of  course,  the  Fox  Island  signal  (in  this  case 
S 2  ( • ) )  is  not  present  and  ( * )  shows  that  some 
activity  is  present  possibly  at  the  second  signal 
velocity. ^However ,  when  the  y^(t)  is  reconstructed, 
say,  by  yt+(t),  the  value  of  the  F  statistic,  which 
must  be  exceeded  in  order  to  reject  the  hypothesis 
that  S2(.)  is  absent  at  the  ,99  probability  level,  is 
2.48.  This  significance  value  is  not  exceeded  at  any 
frequency  as  can  be  seen  from  the  computed  F  statistics 
in  Table  I.  Therefore,  we  accept  the  hypothesis  that 
the  Tonga  signal  is  absent.  Now  a  second  data  array  was 
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F  ’urc*  7.  \na2y*is  of  Tonga  5j(t)  event  using  \  «  19 
c.. .nncls  recorded  at  TFO  at  high  signal  to  noise  ratio. 
Channels  are  nornaliicd  to  equal  peak- to-y rat  aapiitude. 


TABU  1 


Values  of  F  Testing  for  the  Presence  of 
Fox  in  Data  Containing  Only  Tonga 


Ktfcps J 

.156 

.313 

.469 

.625 

.781 

.938 

1.094 

1.250 

1.406 

1.563 

1.719 

1.875 


F 

1.89730 

.99805 

1.96731 

2.10992 

.89259 

.97330 

1.02410 

.98761 

1.02226 

.95154 

.61049 

.60003 


■adc  up  which  contained  the  signals  from  both  Tonga  and 
I  ox  island,  with  the  results  given  in  Figure  8  and  Table 
II.  ke  note  fro*  Figure  8  that  good  reproductions  of  both 
signals  arc  obtained  as  well  as  a  good  reconstructed  tra*e 
y**(t).  The  F  statistic  rejects  the  absence  of  the  Tonga 
signal  at  a  very  high  level  of  significance*  Kc  include 
in  Figure  8  the  estiaates  for  the  Tonga  and  Fox  events 
obtained  by  bcaaforaing  the  original  events  before  they 
arc  aixed.  These  can  serve  as  reference  signals  against 
which  to  judge  the  effectiveness  of  future  experiments 
using  mixed  noisy  data*  As  an  example,  we  consider  the 
noisy  data  in  Figure  9  where  neither  signal  is  visible 
on  the  original  trace  y19(t).  The  filtered  traces 
enhance  the  signal  fairly  well  and  the  F  statistic  is 
still  significant  as  can  be  seen  from  Table  III.  Thus, 
it  is  reasonable  to  assume  that  a  signal  (Fox)  inbeddod 
in  the  code  of  another  signal  (Tonga)  can  be  defected 
using  the  full  array. 

In  order  to  examine  a  case  where  the  effects 
of  distortion  can  be  seen  on  the  straight  beam,  the 
roles  of  the  Tonga  and  Fox  Island  signals  were  inter¬ 
changed  within  the  prograu.  In  Figure  10a 
we  sec  that  the  bean  formed  estimate  S 2 ( • )  contains  a 
definite  contribution  from  the  initial  cycles  of  the 
Tonga  event,  while  the  maximum  likelihood  filters 
reject  this  component.  This  shows  that  the  distortion 
predicted  by  Figure  2  can  be  a  real  factor  for  seismic 
signals.  For  the  purpose  of  future  reference,  we  give 
the  values  of  the  Y  statistic  in  Table  IV. 

The  strong  component  of  the  Tonga  signal  at  1.25  II; 
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recorded  .it  TFO  at 


TABU;  11 


Values  of  F  Testing  for  the  Presence  of  the  Fox  Island 
Signal  in  Data  Containing  Both  Fox  Island  and  Tonga  Signal 


Hz (cps)  F 

.156 

.84619 

.313 

45.39547 

.469 

26.19131 

.625 

58.49830 

.781 

10.86142 

.938 

24.44770 

1.094 

14.68781 

1.2S0 

13.02673 

1.406 

S. 83217 

1.563 

5.08383 

1.719 

2.93364 

1.875 

2.95517 

2,031 

1.54600 

alys  i 


TABU-  III 


Values  of  F  Testing  for  the  Presence  of  the  Fox 
in  Noisy  Data  Containing  Both  Fox  Island  and 


Island  Signal 
Tonga  Signals 


Hz(cps) 

.156 

.313 

.469 

.625 

.781 

.938 

1.094 

1.250 

1*406 

1.563 

1.719 

1.875 

2.031 


F 

6.82573 

2.74730 

2.82802 

14.48280 

6.45614 

8.83391 

5.37031 

3.99074 

1.34506 

1.24787 

1.51457 

.40499 

.91461 


for  N  -  19  channels  showing  bias  in  beam  lormcu  estimate 

for  Fox.  .  ,  ,  . . .  . 

Channels  arc  normalized  to  equal  pcak-to-pcak  amplitude. 


shows  up  in  the  value  of  the  F  statistic  at  1.25  Hz.  An 
even  greater  improvement  over  signal  beamforming  can  be 
seen  in  Figure  10b  in  which  the  data  arc  from  the  seven 
center  elements  at  TFO. 

In  order  to  examine  further  the  detection  and  esti¬ 
mation  capability  of  the  two-signal  model  for  a  some¬ 
what  smaller  array,  we  again  chose  the  seven  center 
elements  at  TFO  and  investigated  the  process  of  testing 
for  Tonga  in  the  presence  of  Fox.  A  number  of  possible 
azimuths  were  tried  for  Tonga  assuming  that  Fox  was  fixed 
at  17  km/scc  and  308  degrees.  This  enables  one  to  examine 
the  sensitivity  of  the  F  statistic  in  resolving  the 
azimuth  of  the  Tonga  signal  in  much  the  same  way  as  a 
frequency  wave  number  spectrum.  Figures  11  and  12  show 
the  high  resolution  frequency  wave  number  spectra  for 
the  seven  center  elements,  for  the  ordinary  and  noisy 
eases  corresponding  to  the  data  in  Figures  8  and  9 
respectively.  Figure  11  shows  that  Tonga  and  Fox  may  be 
resolved  fairly  well  by  conventional  methods  into  two 
separate  components.  Figure  12  shows  that  the  ability 
to  distinguish  two  signals  is  diminished  considerably 
by  noisy  data.  Therefore,  we  might  consider  searching  a 
reasonable  collection  of  azimuths  to  look  at  the  esti¬ 
mation  and  resolution  capabilities  of  the  maximum  likeli¬ 
hood  procedure.  Figure  13  shows  the  estimates  for  the 
waveform  of  Tonga  which  would  be  obtained  for  various 
possible  assumed  azimuths  if  the  seven  center  elements 
at  TFO  were  used  on  data  with  a  high  signal  to 
noise  ratio.  By  comparing  the  waveform  obtained 
with  the  true  version  in  Figure  8^  we  note  that 
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Figure  11.  Frequency  wave  number  analysis  N 
measuring  Tonga  and  Fox  at  TFO  (1.25  Hz). 
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Figure  12.  Frequency  wave  number  analysis  N=  7  noisy 
channels  measuring  Tonga  and  Fox  at  TFO  (1.25  Hz). 
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Figure  13.  Maximum  likelihood  estimate  for  Tonga  using 
Tonga-Fox  mixture  in  N  -  7  channels  at  TFO  (Figure  8) 
with  various  possible  azimuths  for  Tonga. 

Channels  are  normalized  to  equal  peak-to-peak  amplitude. 
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reasonable  estimates  are  obtained  using  any  azimuth  in 
the  range  180-300  degrees.  Definitive  information  on 
the  detector  is  available  in  the  plot  of  the  1  statistic 
against  azimuth  shown  in  Figure  14.  In  this  ease,  values 
of  F  at  the  frequencies  1.25  Hz  and  1.41  Hz  exceed  the 
.001  significance  level  (i.e.  we  would  expect  only  1 
false  alarm  per  1000  frequencies  or  per  1000  reported 
experiments  at  one  frequency)  for  azimuths  between  180 
and  280  degrees  with  the  peak  occurring  between  230  and 
240  degrees.  This  agrees  with  the  approximate  azimuth 
read  from  the  F-K  plot  (Figure  11).  The  two  procedures 
appear  to  work  equally  well  in  this  case  as  the  width 
of  the  main  peak  is  about  10  degrees  in  either  case. 

No  comparisons  will  be  made  on  the  relative 
dropoff  in  db  for  these  procedures  since  the  important 
measures  in  detection  are  not  the  absolute  units  in 
which  a  test  statistic  is  expressed,  but  rather  the 
probability  that  the  test  statistic  would  exceed  the 
speified  threshold  when  the  signal  is  not  there.  For 
example,  the  comparison  of  the  capability  of  a  detector 
which  uses  peak  to  peak  amplitude  with  a  detector  which 
computes  the  ratio  of  the  signal  power  to  the  noise 
power  cannot  be  made  on  the  basis  of  the  values  for 
these  two  test  statistics,  but  must  be  made  on  the 
basis  of  their  probability  distributions  under  presence 
or  absence  of  the  signal. 

A  noisy  case  (Figure  9)  produces  for  the  seven 
channel  ease  the  estimates  for  the  Tonga  waveform  shown 
in  Figure  15.  In  this  case,  the  only  evidence  of  the 
Tonga  signal  is  a  slight  downward  deflection  at  the 
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Figure  IS.  Maximum  likelihood  estimate  for  Tonga  using 
Tonga-Fox  mixture  in  N  -  7  noisy  channels  at  TFO  with 
various  possible  azimuths  for  Tonga. 

Channels  arc  normalized  to  equal  peak-to-pcak  amplitude 


point  where  Tonga  should  appear.  The  plot  of  the  F 
statistic,  Figure  16  however,  clearly  exhibits  a  peak 
(significant  at  .001  level)  at  220  degrees 
azimuth*  In  fact,  the  F  statistic  exceeds  the  .01 
thicshold  between  190  and  270  degrees  as  in  the  high 
signal  to  noise  ratio  ease. 
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AZIMUTH  (DEGREES) 


Figure  16.  F  statistic  for  detecting  Tonga  in  Tonga-Fox 
noisy  mixture  (N  ■  7  center  elements  at  TFO). 


SUMMARY 


An  investigation  has  been  made  of  the  effec¬ 
tiveness  of  truncated  likelihood  filters  in  simul¬ 
taneously  estimating  and  detecting  two  real  seismic 
signals  mixed  with  noise.  It  is  found  that  the  maximum 
likelihood  filters  give  better  estimates  than  simple 
beam  forming,  with  the  superiority  most  pronounced  in 
small  arrays.  A  theoretical  procedure  for  examining 
the  distortion  of  maximum  likelihood  and  bcamformed 
filters  is  derived  and  illustrated.  The  detection 
capabilities  of  the  F  statistic  are  examined  for  a  small 
seven  clement  short  period  array  at  TFO.  A  comparison 
with  high  resolution  FK  spectra  shows  that  the  two- 
signal  version  of  the  F  statistic  may  be  superior  in 
the  case  where  the  noise  level  is  high. 
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